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Two  algorithms  are  described  for  assessing  systematic  errors  in  acoustic  travel-time  tomography  of 
the  atmosphere,  the  goal  of  which  is  to  reconstruct  the  temperature  and  wind  velocity  fields  given 
the  transducers’  locations  and  the  measured  travel  times  of  sound  propagating  between  each 
speaker-microphone  pair.  The  first  algorithm  aims  at  assessing  the  errors  simultaneously  with  the 
mean  field  reconstruction.  The  second  algorithm  uses  the  results  of  the  hist  algorithm  to  identify 
the  ray  paths  corrupted  by  the  systematic  errors  and  then  estimates  these  errors  more  accurately. 

Numerical  simulations  show  that  the  first  algorithm  can  improve  the  reconstruction  when  relatively 
small  systematic  errors  are  present  in  all  paths.  The  second  algorithm  significantly  improves  the 
reconstruction  when  systematic  errors  are  present  in  a  few,  but  not  all,  ray  paths.  The  developed 
algorithms  were  applied  to  experimental  data  obtained  at  the  Boulder  Atmospheric  Observatory. 
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I.  INTRODUCTION 

The  goal  of  acoustic  travel-time  tomography  is  to  recon¬ 
struct  physical  properties  of  the  medium  given  the  coordi¬ 
nates  of  speakers  and  microphones  and  the  travel  times  of 
sound  impulses  propagating  through  the  medium.  In  tomog¬ 
raphy  of  the  atmosphere,  the  medium  properties  of  interest 
are  the  temperature  and  wind  velocity  fields  within  an  area 
(two-dimensional  problem) 1-6  or  volume  (three-dimensional 
problem).7'8  Mathematically,  this  problem  belongs  to  the 
class  of  ill-posed  problems;  namely,  the  solution  may  not  be 
unique  for  a  given  set  of  data.  For  example,  some  wind  veloc¬ 
ity  fields  with  non-zero  divergence  (i.e.,  with  volume  sources 
or  sinks)  may  not  affect  travel  times  along  the  sound  propa¬ 
gation  paths.9-12  Moreover,  there  may  exist  temperature  and 
wind  velocity  fluctuations  that  are  too  small  to  be  resolved, 
given  the  spatial  sampling  of  the  available  paths.  Such  fields 
are  undetectable  by  travel-time  tomography  without  supple¬ 
mentary  data  and  are  not  considered  in  this  paper. 

The  need  to  properly  estimate  and  eliminate  systematic 
errors  is  encountered  in  many  disciplines.  Estimation  algo¬ 
rithms  can  be  grouped  into  two  categories:  (1)  algorithms  for 
the  estimation  and  prediction  of  the  systematic  errors  owing 
to  incompleteness  of  the  numerical  model  of  a  physical  phe¬ 
nomena  being  modeled13-16  and  (2)  statistical  algorithms  for 
the  systematic  errors  in  the  measured  data.1718  The  focus  of 
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this  paper  is  on  the  second  category,  when  reconstruction  of 
the  atmospheric  fields  is  done  using  travel-time  tomography. 

In  acoustic  tomography  experiments,  the  transducers’ 
coordinates  and  the  travel  times  of  sound  propagation 
between  speakers  and  microphones  should  be  measured 
accurately.  The  travel-time  measurements  include  assess¬ 
ment  of  the  time  delays  in  electronic  circuits  and  mechanical 
hardware  (e.g.,  drivers  and  microphones)  of  a  tomography 
array.  References  4  and  5  already  discussed  the  effect  of  ran¬ 
dom  errors  in  the  travel  times  on  the  tomographic  recon¬ 
struction.  The  present  paper  considers  systematic  errors  in 
the  travel  times  that  may  be  significantly  larger  than  random 
errors  and,  thus,  can  significantly  affect  the  quality  of  the 
reconstruction  of  temperature  and  wind  velocity  fields. 

Several  factors  can  contribute  to  the  systematic  errors  in 
the  travel-time  measurements  of  acoustic  tomography  of  the 
atmosphere.  First,  the  systematic  errors  are  caused  by  the 
errors  in  measurements  of  the  time  delays  of  signal  propaga¬ 
tion  in  hardware  and  electronic  circuits  of  the  tomography 
array  and  errors  in  synchronization  of  the  transmitted  and 
recorded  signals.  For  example,  if  synchronization  is  done 
using  the  computer  operating  system,  the  corresponding  sys¬ 
tematic  error  can  be  of  order  1  ms.  Second,  the  errors  in 
transducers  coordinates  can  be  as  large  as  30  cm.  These 
errors  are  equivalent  to  the  systematic  errors  in  the  travel 
times  of  0.9  ms.  Third,  loudspeakers  which  are  used  in 
acoustic  tomography  are  extended  sources.  To  simplify  a  to¬ 
mography  experiment  with  many  speakers  and  microphones 
(or  if  a  tomography  array  needs  to  be  assembled  quickly), 
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this  fact  can  be  ignored  in  the  mathematical  formulation  and 
loudspeakers  can  be  considered  as  point  sources.  This 
approximation  can  result  in  the  systematic  errors  of  the  coor¬ 
dinates  of  the  effective  point  sources  of  order  1  m,  which  is 
equivalent  to  3  ms.  Finally,  there  might  be  other  sources  of 
systematic  errors  which  are  difficult  to  identify  but  which 
can  be  as  large  as  those  mentioned  above. 

The  systematic  errors  might  slowly  change  in  time 
because  of  changing  conditions  in  electronic  circuits,  me¬ 
chanical  hardware,  and  the  positions  of  transducers.  The  im¬ 
portance  of  error  analysis  in  tomography  problems  is  well 
understood  and  the  assessments  of  their  influence  on  the 
reconstructions  are  available.  ’  ’  The  current  paper  is 
focused  on  the  techniques  to  mitigate  these  errors,  which  de¬ 
spite  all  efforts  for  their  elimination  are  still  present  in  the 
data.  The  paper  describes  two  algorithms  for  an  explicit 
assessment  of  the  systematic  errors  from  the  corrupted 
travel-time  measurements  and  discusses  the  effect  of  these 
systematic  errors  on  the  quality  of  the  reconstruction.  First, 
the  algorithms  are  tested  on  numerically  simulated  tempera¬ 
ture  and  wind  velocity  fields  using  large  eddy  simulation 
(LES).22  Then,  by  applying  the  algorithms  to  experimental 
data  obtained  at  the  Boulder  Atmospheric  Observatory 
(BAO),23  the  paper  shows  how  the  algorithms  can  improve 
tomographic  reconstructions. 


II.  PROBLEM  STATEMENT 

The  basic  equation  for  the  reconstruction  of  temperature 
and  wind  velocity  fields  can  be  found,  for  example,  in  Refs. 
4  and  7.  In  particular,  the  two-dimensional,  linearized  prob¬ 
lem  is  described  by  the  following  equation: 

uo(t)cos  (pj  +  ro(0sin  Vi 

.  cow  ; 

-F(t)  +  ei(t),  (1) 


where  t{[  is  the  travel  time  along  the  /th  sound  propagation 
path;  /  =  1,  2,...,  /;  t  is  the  time  at  which  the  measurement 
was  performed;  L,  is  the  length  of  the  /th  path;  Cq  is  the  spa¬ 
tial  mean  sound  speed  in  the  motionless  medium;  u0  and  Vq 
are  the  x  and  y  Cartesian  components  of  the  spatial  mean 
two-dimensional  wind  velocity;  <p,-  is  the  angle  between  the 
positive  direction  of  the  x  axis  and  the  direction  of  the  /th 
path;  and  i:,  represents  the  random  errors  in  measurements  of 
the  travel  times.  In  Eq.  (1),  F  accounts  for  travel-time  contri¬ 
butions  caused  by  fluctuations  of  the  temperature  T( r,  t)  and 
wind  velocity  [h( r,  t),  u(r,  t)\  in  the  tomographic  area  rela¬ 
tive  to  T0  and  c0,  respectively, 


F(t)  = 


cl(t) 


dl 


co  (t) 


27b(0 


T(r,  t) 


+«( r,  /)cos  cpi  +  v(r,  /)sin  cp: 


(2) 


where  r  =  (x,  y )  is  the  two-dimensional  space  vector,  the 
integration  is  implemented  along  the  ray  path,  and  T0  is  the 
spatial  mean  temperature  related  to  Cq  by  a  formula 


Cq  =  i'R„To  with  y  «  1 .40  being  the  ratio  of  the  specific  heats 
and  Ra  «  287.058  J/(kg  K)  being  the  specific  gas  constant 
for  dry  air.  The  formulation  given  by  Eqs.  (1)  and  (2) 
assumes  straight-line  sound  propagation,  which  is  valid  if 
the  vertical  gradients  of  temperature  and  wind  magnitude  are 
small,  so  that  the  rays  belong  to  a  single  plane,  and  the  total 
wind  fields  and  temperature  fluctuations  much  smaller  than 
the  mean  sound  speed  and  mean  temperature,  respectively: 
|m0  +  m|/c0<C  1,  |o0+  r|/c0<C  1,  and  \T\/T0<g.l.  Keeping 
only  the  first-order  powers  of  these  small  ratios  and  using 
the  results  in  Sec.  3.3.1  of  Ref.  24,  it  can  be  shown  that 
n  ■  s  =  1.  Therefore,  in  this  approximation,  the  unit  vectors  n 
(normal  to  a  wave  front)  and  s  (direction  of  impulse  propaga¬ 
tion)  coincide.  As  shown  in  Ref.  1,  the  linear  approximation 
in  acoustic  tomography  of  the  atmosphere  introduces  negli¬ 
gible  errors  for  relatively  short  sound  propagation  paths  on 
the  order  of  200  m.  The  temperature  in  Eq.  (2)  represents  the 
acoustic  virtual  temperature,  Tav  =  7),  +  T.  A  thermody¬ 
namic  temperature  Tth  can  be  recalculated  from  the  acoustic 
virtual  temperature  using  the  following  relationship  [Eq. 
(6.23)  in  Ref.  24]:  Tav  =  7)/,(l  +  0.51  \q),  where  q  is  the  spe¬ 
cific  humidity  of  air.  As  discussed  in  Refs.  9-12,  for  a 
unique  reconstruction  of  u{ r,  t),  and  u(r,  t),  we  shall  assume 
that  their  sources  and  sinks  are  negligible  at  each  time  t,  that 
is  du(r1  t)/dx  +  <9r(r,  t)/dy  «  0.  The  actual  assumption 
employed  in  LES  is  du(r,t)  jdx  +  dv (r ,  t)/dy  +  dw(r ,  t)  /dz 
=  0,  as  valid  for  an  incompressible,  three-dimensional  fluid. 
The  analysis  of  the  LES  fields  used  in  this  paper  corrobo¬ 
rated  validity  of  the  small  two-dimensional  divergence 
assumption. 

The  reconstruction  of  the  temperature  and  wind  velocity 
fields  can  be  implemented  in  two  stages,  as  proposed  in 
Ref.  4.  First,  the  mean  values  of  the  sound  speed  and  wind 
velocity  components  are  reconstructed  via  the  least-squares 
fit  using  Eq.  (1)  with  the  F(t)  term  omitted.  This  reconstruc¬ 
tion  will  be  called  “Algorithm  1”  in  this  paper.  Estimation  of 
the  variance  of  the  reconstruction  errors  due  to  measurement 
errors  £,  is  implemented  by  a  standard  error  estimation  tech¬ 
nique25  that  assumes  independence  of  the  measurement 
errors,  both  in  repeated  measurements  for  the  same  sound 
path  and  between  any  two  paths.  Then,  the  F(t)  term  is  cal¬ 
culated  using  Eq.  (1),  which  enters  as  the  input  data  in  the 
left-hand  side  of  Eq.  (2)  for  the  reconstruction  of  the  temper¬ 
ature  and  wind  velocity  fluctuations  using  the  time- 
dependent  stochastic  inversion  (TDSI)  algorithm. 4-7  The 
advantage  of  TDSI,  compared  to  other  inverse  techniques,  is 
that  it  uses  data  measured  multiple  times  to  reconstruct  the 
fluctuations  at  one  particular  time  (effectively  increasing  the 
amount  of  data),  incorporating  known  spatial-temporal  cor¬ 
relations  of  the  fluctuations  into  the  reconstruction. 

The  problem  with  Algorithm  1  is  that  the  travel  times  t’[ 
can  be  corrupted  by  the  systematic  errors  d,  which  are  con¬ 
stant  or  slowly  changing  in  time  and  might  be  different  for 
each  path  /.  There  are  two  principal  distinctions  between  sys¬ 
tematic  and  random  errors.  First,  systematic  errors  are  con¬ 
stant  (or  vary  slowly)  from  one  measurement  to  another 
while  random  errors  are  independent  and  have  zero  mathe¬ 
matical  expectation.  Second,  their  magnitudes  might  be 
larger  than  that  of  £,.  As  demonstrated  in  Sec.  IV,  large 
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systematic  errors  in  the  travel  times  make  the  reconstruction 
of  temperature  and  wind  velocity  impossible.  The  goal  of 
this  paper  is  to  investigate  whether  these  systematic  errors  in 
the  travel  times  can  be  explicitly  estimated  from  the  meas¬ 
urements  and,  consequently,  eliminated  from  the  data. 

III.  METHOD 

Our  procedure  for  estimating  systematic  errors  assumes 
that  these  errors  are  constant  in  a  series  of  repeated  measure¬ 
ments  of  the  travel  times  during  some  time  interval  At.  It  is 
assumed  that  systematic  errors  with  a  weak  trend  can  also  be 
treated  this  way,  if  they  are  nearly  constant  during  the  inter¬ 
val  of  data  measurements  At,  which  was  equal  to  10  s  in  the 
outdoor  and  numerical  experiments.  Let  d,  be  a  systematic 
error  in  the  ith  ray.  Then,  Eq.  (1)  can  be  recast  as 

jr,  a  _  U  A  _  u0  cos  (p;  +  v0  sin  <?A 

'  c0  V  c0  ) 

—  F(t)  +  £;(t)  +  di ,  t  £  At.  (3) 

Unlike  Eq.  (1),  the  mean  values  Co,  Up,  and  Vp,  as  well  as  the 
errors  dh  are  considered  as  constants  during  the  time  interval 


At  for  which  Eq.  (3)  is  applied.  This  does  not  introduce  any 
additional  errors  or  inconsistencies  in  the  reconstruction  algo¬ 
rithm  outlined  above.  The  only  difference  is  that  the  mean 
values  Co,  u0,  and  Vp  are  now  not  instantaneous  spatial  mean 
values  but  the  spatial-temporal  mean  values  over  the  tomo¬ 
graphic  area  defined  by  the  tomographic  setup  and  over  the 
time  interval  At.  This,  of  course,  changes  the  estimation  of 
F(t)  and,  consequently,  the  reconstructed  temperature  and 
wind  velocity  fluctuations.  As  before,  the  reconstruction  of 
these  mean  values  is  implemented  by  neglecting  the  fluctua¬ 
tion  term  F(t)  in  Eq.  (3).  However,  unlike  the  original 
Algorithm  1,  the  mean  values  are  estimated  simultaneously 
with  the  unknown  systematic  errors  dh  using  the  least- 
squares  fit.  Such  an  idea  was  exploited  in  underwater  acoustic 
tomography  in  a  variety  of  reconstruction  techniques. 26-28 
This  reconstruction  will  be  referred  to  as  “Algorithm  2.” 
Note  that  one  needs  to  use  data  from  several  measurements 
to  make  the  problem  over  determined.  Indeed,  a  set  of  travel¬ 
time  measurements  implemented  at  some  time  t\  €  At  has  7 
data  and  7+3  unknowns  (co,  u0,  v0,  and  all  d,’s).  Using  data 
from  another  consecutive  measurement  at  U  €  At,  one  has  27 
data  and  still  7+3  unknowns.  Thus,  for  7  >  3,  one  should 
use  at  least  two  data  sets  acquired  at  two  different  times.  To 
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FIG.  1.  (Color  online)  LES  fields  used  in  the  numerical  study,  (a)  Temperature  field.  Black  streamlines  with  arrows  indicate  wind  direction,  (b)  First  horizontal 
component  of  the  wind  velocity  vector,  u0  +  u.  (c)  Second  horizontal  component  of  the  wind  velocity  vector,  v0  +  v.  (d)  Principal  schema  of  the  tomography 
sensor  array. 
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(a)  (b) 


Sound  path  index  Sound  path  index 

FIG.  2.  (Color  online)  (a)  Travel  times  for  the  15  paths  shown  in  Fig.  1(d).  (b)  Actual  (asterisks)  and  estimated  (circles)  systematic  errors  in  the  travel  times. 
Note  that  the  actual  systematic  errors  are  zeros  for  this  case. 


increase  the  robustness  and  reduce  the  errors  of  the  estima¬ 
tion,  more  data  can  be  used.  The  reconstmction  of  the  fluctu¬ 
ations  is  then  performed  by  the  TDSI  algorithm,  as  before.  In 
the  next  section,  the  proposed  method  is  studied  in  numerical 
experiments  with  simulated  temperature  and  wind  velocity 
fields. 

IV.  NUMERICAL  STUDY 

For  this  study,  atmospheric  temperature  and  wind  veloc¬ 
ity  fields  were  created  using  LES22  that  represents  realistic 
atmospheric  fields.  In  LES,  the  Navier-Stokes  equations  are 
solved  numerically  for  relatively  large  spatial  scales  (in  our 
case,  larger  than  4  m),  while  the  fields  at  smaller  scales  are 
parameterized  depending  on  the  regime  of  the  atmosphere 
and  ground  conditions.  In  the  present  study,  a  two- 
dimensional  slice  (at  a  height  of  13.75  m)  of  the  fields  corre¬ 
sponding  to  a  typical  unstable  atmosphere  over  a  homogene¬ 
ous  ground  was  used.  The  frozen  turbulence  hypothesis29,30 
was  used  to  create  temporal  evolution  of  turbulence  from  the 
two-dimensional  slice  of  the  LES  fields.  According  to  this 
hypothesis,  turbulence  translates,  without  distortion,  at  a  con¬ 
stant  advection  velocity.  In  our  numerical  simulations,  the 
advection  velocity  was  set  to  have  Cartesian  components  of 
(4,  4)  m/s.  These  time-dependent  fields  allowed  us  to  calcu¬ 
late  travel  times  t‘[  for  different  times  t  simulating  physical 
measurements.  A  realization  of  the  temperature  and  wind  ve¬ 
locity  fields  to  be  reconstructed  and  the  tomographic  setup 
are  shown  in  Fig.  1.  Black  streamlines  with  arrows  in  Fig. 
1(a)  indicate  the  direction  of  wind.  Figures  1(b)  and  1(c) 
depict  the  components  of  the  wind  velocity,  u0  +  u  and  i;0  + 
v,  respectively.  Figure  1(d)  shows  the  tomographic  setup  con¬ 
sisting  of  five  microphones  and  three  speakers.  This  figure  is 
a  view  from  above  on  the  BAO  acoustic  tomography  array. 
For  a  single  scan,  the  setup  provides  15  travel  times  corre¬ 
sponding  to  15  rays.  The  travel  times  t'[  corresponding  to  the 
fields  in  Fig.  1  are  shown  in  Fig.  2(a)  versus  the  index  i  of  a 
particular  path.  In  the  numerical  simulations,  the  contribu¬ 
tions  to  the  travel  times  due  to  fluctuations  [i.e.,  the  term  F  in 
Eqs.  ( 1 ) — (3)]  were  on  the  order  of  ±4ms.  To  be  able  to 
reconstmct  the  fluctuations,  we  need  the  accuracy  of  their 


measurements  to  be  less  than  10%.  That  is,  a  threshold  of 
0.4  ms  was  imposed  as  the  acceptable  error  magnitude  in 
travel  times.  Since  the  experimental  setup  at  BAO  was  sup¬ 
posed  to  provide  this  accuracy,  this  value  is  also  a  threshold 
for  deciding  whether  systematic  errors  are  present  in  the  data. 
The  following  cases  are  considered  to  investigate  whether 
large  but  constant  systematic  errors  can  be  effectively  esti¬ 
mated  and  eliminated  from  the  data. 

A.  Case  1 

There  are  no  systematic  errors  in  the  data.  This  case 
illustrates  the  situation  when  the  data  processing  algorithm 
assumes  the  presence  of  the  systematic  errors,  the  least- 
squares  fit  is  based  on  Eq.  (3),  but  there  are  no  systematic 
errors  in  reality.  As  shown  in  Fig.  2(b),  the  least-squares  fit 
(Algorithm  2)  yields  artifact  systematic  errors  (circles)  rang¬ 
ing  from  —0.20  ms  to  0.23  ms.  For  this  estimate,  travel  times 
from  five  consecutive  scans  were  used,  providing  75  data 
points.  As  one  can  see  in  Table  I,  despite  spurious  recon¬ 
structed  systematic  errors,  the  reconstruction  of  the  mean 
fields  is  accurate,  as  good  as  by  the  original  Algorithm  1 
designed  for  no  systematic  error  reconstruction.  Both  algo¬ 
rithms  yield  close  values  and  comparable  estimated  errors  of 
the  reconstruction. 

B.  Case  2 

There  is  a  systematic  error  in  a  single  ray.  For  this  case, 
a  constant  systematic  error  of  10  ms  was  added  to  the  travel 
time  of  the  sound  impulse  propagating  along  the  third  path, 
f  =  3,  which  introduces  a  5%  error  in  the  travel  time.  Note 

TABLE  I.  Mean  field  reconstructions  without  systematic  errors  in  the  input 
data.  “Algorithm  1”  and  “Algorithm  2”  refer  to  the  original  algorithm  with¬ 
out  systematic  errors  and  the  proposed  one,  assuming  systematic  errors  in 
all  paths,  respectively. 


Mean  fields 

T(°C) 

u  (m/s) 

v  (m/s) 

Actual 

28.53 

4.00 

4.00 

Algorithm  1 

28.60  ±  0.09 

4.03  ±  0.07 

4.05  ±  0.07 

Algorithm  2 

28.50  ±  0.05 

4.03  ±  0.04 

4.10  ±0.04 
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FIG.  3.  (Color  online)  Actual  (asterisks)  and  reconstructed  systematic  errors.  Circles  depict  the  Algorithm  2  estimate,  assuming  that  systematic  errors  are  pres¬ 
ent  in  all  rays.  Squares  depict  the  estimation  by  Algorithm  3  when  the  path  indices  with  systematic  errors  are  determined  from  Algorithm  2.  (a)  One  systematic 
error  in  path  3.  (b)  Two  systematic  errors  in  paths  3  and  9. 


that  10  ms  is  of  order  of  the  combined  systematic  errors  men¬ 
tioned  in  Sec.  I.  Such  a  systematic  error  is  prohibitively  large 
for  reconstruction  of  weak  fluctuations.  Without  proper 
detection  and  elimination  of  such  errors,  the  reconstruction 
will  be  severely  distorted,  as  demonstrated  below.  In  the  out¬ 
door  experiments.  Sec.  V,  the  estimated  systematic  errors 
are  on  the  order  of  1  ms.  Figure  3(a)  depicts  the  actual  sys¬ 
tematic  errors  (asterisks)  and  the  estimated  systematic  errors 
obtained  with  Algorithm  2  (circles),  assuming  that  the  errors 
exist  for  each  ray.  As  one  can  see  from  Fig.  3(a),  this  algo¬ 
rithm  does  not  provide  a  reliable  estimate.  Namely,  the  sys¬ 
tematic  error  in  a  single  ray  causes  spurious  reconstructed 
systematic  errors  in  all  rays.  Although  the  third  (true)  error 
is  the  largest,  7.2  ms,  this  is  not  a  good  estimation  of  the 
actual  10  ms  error;  their  discrepancy  is  significantly  larger 
than  the  imposed  threshold  of  0.4  ms.  Additionally,  the  other 
errors  are  not  of  negligible  magnitudes,  ranging  from 
—3.2  ms  to  1.0  ms.  Moreover,  the  mean  field  reconstructions 
by  both  algorithms  become  inaccurate;  see  “Algorithm  1” 
and  “Algorithm  2”  in  Table  II.  The  discrepancy  between  the 
actual  and  reconstructed  mean  fields  is  as  large  as  2.7  °C, 
1.1  m/s,  and  2.6m/s  for  Algorithm  1  and  3.6°C,  1.2 m/s,  and 
3.0  m/s  for  Algorithm  2.  Note  that  the  estimated  reconstruc¬ 
tion  errors  in  the  Algorithm  1  have  increased  by  1  order  of 
magnitude  in  comparison  with  case  1  when  the  travel  times 
were  error  free.  Therefore,  such  large  estimated  errors  can 
serve  as  indicators  of  the  presence  of  systematic  errors  in  the 
data. 

TABLE  II.  Mean  field  reconstructions  with  a  single  systematic  error. 
“Algorithm  1”  and  “Algorithm  2”  refer  to  the  algorithms,  without  systematic 
errors  and  assuming  systematic  errors  in  all  paths,  respectively.  "Algorithm 
3”  refers  to  the  refined  estimation  algorithm,  assuming  a  systematic  error  in 
ray  3. 


Mean  fields 

T(°C) 

u  (m/s) 

v  (m/s) 

Actual 

28.53 

4.00 

4.00 

Algorithm  1 

25.9  ±  0.9 

5.1  ±0.7 

6.6  ±0.7 

Algorithm  2 

24.95  ±  0.04 

5.22  ±  0.04 

7.03  ±  0.03 

Algorithm  3 

28.67  ±  0.09 

4.01  ±  0.07 

3.99  ±  0.08 

1806  J.  Acoust.  Soc.  Am.,  Vol.  134,  No.  3,  September  2013 


To  mitigate  these  problems,  we  propose  a  refined  recon¬ 
struction,  referred  to  as  “Algorithm  3”  hereafter.  The  idea  is 
to  use  the  results  of  the  Algorithm  2  to  determine  which  paths 
actually  have  the  errors.  Judging  from  Fig.  3(a),  the  largest 
estimated  error  occurs  in  the  third  path.  Now,  the  least- 
squares  estimation  can  be  repeated  analogously  to  Algorithm 
2  with  systematic  errors,  appearing  in  Eq.  (3),  but  only  a  sin¬ 
gle  systematic  error  is  accounted  for  in  the  third  path.  The 
systematic  errors  of  a  reconstruction  refined  in  this  manner 
are  shown  in  Fig.  3(a)  by  squares,  and  the  reconstructed  mean 
fields  are  shown  in  Table  II.  The  refined  reconstruction  yields 
significantly  better  results.  First,  the  estimated  error  in  the 
third  ray  is  10.2  ms,  which  is  very  close  to  the  actual  value  of 
10  ms,  with  the  discrepancy  below  the  0.4  ms  threshold. 
Second,  the  mean  field  values  are  significantly  better  than  in 
the  two  previous  algorithms.  The  differences  between  the 
actual  and  reconstructed  fields  are  small:  0.14°C,  0.1  m/s,  and 
0.1  m/s.  The  estimated  reconstruction  errors  have  also 
decreased  as  compared  to  the  original  reconstruction  by 
Algorithm  1.  To  validate  the  refined  reconstruction  approach 
(Algorithm  3),  let  us  consider  the  following  case. 

C.  Case  3 

There  are  distinct  systematic  errors  in  two  rays.  The  sys¬ 
tematic  errors  of  10  and  15  ms  were  added  to  the  third  and 
ninth  rays,  respectively,  which  corresponds  to  the  5.0%  and 
6.5%  relative  errors  in  the  third  and  ninth  travel  times.  Figure 
3(b)  shows  the  actual  (asterisks)  and  estimated  (circles  and 

TABLE  III.  Mean  field  reconstructions  with  systematic  errors  in  two  paths. 
“Algorithm  1”  and  “Algorithm  2”  refer  to  the  algorithms  without  systematic 
errors  and  assuming  systematic  errors  in  all  paths,  respectively.  “Algorithm 
3”  refers  to  the  refined  estimation  algorithm  assuming  systematic  errors  in 
two  paths. 


Mean  fields 

T(°C) 

u  (m/s) 

v  (m/s) 

Actual 

28.53 

4.00 

4.00 

Algorithm  1 

24  ±  1 

3  ±  1 

5  ±  1 

Algorithm  2 

21.09  ±0.04 

0.96  ±  0.04 

5.24  ±  0.03 

Algorithm  3 

28.66  ±0.10 

4.01  ±  0.07 

3.99  ±  0.08 
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squares)  systematic  errors.  As  in  Fig.  3(a),  the  circles  stand 
for  Algorithm  2  (the  two  largest  estimates  indicating  the  path 
indices  with  the  systematic  errors);  and  the  squares  depict  the 
refined  reconstruction  of  the  systematic  errors.  Similarly  to 
case  2,  the  refined  reconstmction  yields  very  accurate 


estimations:  10.2  ms  for  the  actual  10  ms  in  the  third  travel 
time,  and  15.0  ms  for  the  actual  15  ms  in  the  ninth  travel  time. 
The  reconstruction  of  the  mean  fields  by  the  three  algorithms 
is  summarized  in  Table  III.  With  the  systematic  errors  in  two 
data  points.  Algorithm  1  and  Algorithm  2  yield  spurious 


FIG.  4.  (Color  online)  Reconstmction  of  the  unstable  LES  fields  in  the  presence  of  two  constant  systematic  errors  in  paths  3  and  9.  Note  significant  improve¬ 
ment  in  reconstructions  (right  column)  when  Algorithm  3  was  used,  correctly  assuming  systematic  errors  in  paths  3  and  9.  (a)  Temperature  field,  no  estimation 
of  the  systematic  errors,  (b)  Temperature  field,  refined  estimation  of  the  systematic  errors,  (c)  First  horizontal  component  of  the  wind  velocity  vector,  u0  +  u, 
no  estimation  of  the  systematic  errors,  (d)  First  horizontal  component  of  the  wind  velocity  vector,  Uq  +  u,  refined  estimation  of  the  systematic  errors,  (e) 
Second  horizontal  component  of  the  wind  velocity  vector,  v0  +  v,  no  estimation  of  the  systematic  errors,  (f)  Second  horizontal  component  of  the  wind  velocity 
vector,  Vq  +  v,  refined  estimation  of  the  systematic  errors. 
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FIG.  5.  (Color  online)  Actual  (asterisks)  and  estimated  by  Algorithm  2  (circles)  systematic  errors  present  in  all  paths,  (a)  Large  systematic  errors  in  the  interval 
[—5  5]  ms.  (b)  Small  systematic  errors  in  the  interval  [—1  1]  ms. 


estimates.  In  contrast.  Algorithm  3  yields  not  only  the  correct 
systematic  errors  but  also  accurate  reconstructions  of  the 
mean  fields  (see  Table  III).  Moreover,  the  reconstruction  of 
the  fluctuations  and,  thus,  complete  fields  (Fig.  4)  from  these 
distorted  data  using  the  refined  mean  field  reconstruction  is 
also  very  good  despite  mismatched  eddy  distributions  of  the 
temperature  field:  Comparison  of  Fig.  1  and  Fig.  4  reveals 
that  Algorithm  3  yields  reconstructions  with  average  errors  of 
order  0.5  °C  and  0.5  m/s  while  Algorithm  1  yields  spurious 
reconstructions  with  an  error  order  of  7  °C  and  7  m/s.  Figures 
4(a),  4(c),  and  4(e)  illustrate  the  spurious  solution  obtained  by 
Algorithm  1  followed  by  TDSI  for  the  reconstruction  of  the 
fluctuations  from  the  two-path-error  data.  All  fields  are  out  of 
plausible  ranges,  and  the  shape  of  the  fluctuations  mismatches 
those  of  the  actual  fields  shown  in  Fig.  1.  In  contrast, 
Algorithm  3  followed  by  TDSI,  Figs.  4(b),  4(d),  and  4(f), 
yields  much  better  reconstructions,  not  only  in  the  right  range 
of  values  but  also  by  capturing  the  actual  shapes  of  some  fluc¬ 
tuations;  see,  for  example,  the  reconstruction  of  the  u 
component. 

D.  Case  4 

There  are  systematic  errors  in  all  rays.  For  this  case,  sys¬ 
tematic  errors  randomly  distributed  in  the  range  [—5,  5]  ms 
were  added  to  all  travel  times,  introducing  relative  errors  in 
the  travel  times  ranging  from  0.2%  to  4%.  Note  that 
Algorithm  3,  in  this  case,  coincides  with  Algorithm  2.  The 
actual  and  reconstructed  systematic  errors  are  shown  in  Fig. 
5(a).  None  of  the  errors  are  reconstructed  accurately  enough 
resulting  in  a  spurious  reconstruction  of  the  mean  fields;  see 
Table  IV,  rows  “Alg.  1,  large  errors”  and  “Alg.  2,  large 
errors.”  These  results  can  be  summarized  as  follows:  If  all 
input  data  to  a  reconstruction  algorithm  have  large  unknown 
systematic  errors,  as  in  the  considered  case  4,  the  correct 
reconstruction  is  impossible.  However,  if  systematic  errors 
are  of  small  magnitude,  this  somewhat  improves  the  situa¬ 
tion.  Figure  5(b)  shows  actual  and  reconstructed  systematic 
errors  when  the  actual  ones  are  in  the  range  [  —  1,  1]  ms  (rela¬ 
tive  errors  in  the  travel  times  do  not  exceed  1.2%).  In  this 
case,  the  reconstruction  of  the  mean  fields  becomes  better,  as 
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one  can  see  in  Table  IV,  “Alg.  1,  small  errors”  and  “Alg.  2, 
small  errors”  rows.  Note  that  Algorithm  2  yields  slightly  bet¬ 
ter  estimations  than  the  original  Algorithm  1 . 

From  the  considered  cases,  we  conclude  that  system¬ 
atic  errors  in  several,  but  not  all,  travel  times  can  be  suc¬ 
cessfully  detected  and  eliminated  by  the  refined 
reconstruction  technique.  Algorithm  3.  Large  estimated 
errors  in  the  mean  field  reconstructions  obtained  by  the 
original  Algorithm  1  indicate  possible  presence  of  system¬ 
atic  errors  in  the  data.  If  systematic  errors  are  present  in  all 
data,  their  accurate  estimation  and  mean  field  reconstruc¬ 
tion  are,  generally,  impossible.  However,  for  small  system¬ 
atic  errors,  Algorithm  2  may  provide  better  reconstruction 
of  the  mean  fields  than  the  original  Algorithm  1 . 

Algorithms  1 ,  2,  and  3  were  also  used  in  numerical  anal¬ 
ysis  of  the  systematic  errors  for  other  LES  fields  correspond¬ 
ing  to  different  snapshots  of  an  unstable  atmosphere,  a  very 
unstable  atmosphere,  and  a  neutral  atmosphere.  The  results 
supported  the  tendencies  described  above.  For  example. 
Table  V  shows  the  mean  field  reconstructions  for  all  cases 
studied  above  for  very  unstable  LES,  depicted  in  Figs.  6(a), 
6(b),  and  6(c).  One  can  observe  the  same  features  as  for  the 
unstable  LES:  For  errors  in  a  single  path  and  two  paths, 
Algorithm  3  outperforms  both  Algorithm  1  and  Algorithm  2. 
If  large  (±5  ms)  errors  are  present  in  all  paths,  then  all  algo¬ 
rithms  fail  to  reconstruct  the  fields.  However,  if  the  errors  are 
small  (±lms),  the  reconstruction  converges  to  the  actual 
fields,  and  Algorithm  2  slightly  outperfonns  Algorithm  1. 

These  features  are  also  clearly  seen  in  Fig.  6.  Figures 
6(d),  6(e),  and  6(f)  depict  failed  reconstructions  of  temperature 


TABLE  IV.  Mean  field  reconstructions  with  systematic  errors  in  all  input 
data.  “Alg.  1”  and  “Alg.  2”  refer  to  the  algorithms  without  systematic  errors 
and  assuming  systematic  errors  in  all  paths,  respectively. 


Mean  fields 

T  (°C) 

u  (m/s) 

v  (m/s) 

Actual 

28.53 

4.00 

4.00 

Alg.  1 ,  large  errors 

22  ±  1 

3.7  ±0.9 

5  ±  1 

Alg.  2,  large  errors 

23.56  ±  0.04 

4.53  ±  0.04 

5.74  ±  0.03 

Alg.  1 ,  small  errors 

27.8  ±  0.2 

4.3  ±0.2 

4.8  ±0.2 

Alg.  2,  small  errors 

28.02  ±  0.05 

3.94  ±  0.04 

4.35  ±  0.04 
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FIG.  6.  (Color  online)  Reconstruction  of  a  very  unstable  LES.  (a)  Actual  temperature  field.  Black  streamlines  with  arrows  indicate  wind  direction,  (b)  Actual 
first  horizontal  component  of  the  wind  velocity  vector,  Uq  +  u.  (c)  Actual  second  horizontal  component  of  the  wind  velocity  vector,  Vo  +  v.  Reconstruction 
without  estimation  of  the  systematic  errors  (Algorithm  1)  with  two  travel  times  present  in  paths  3  and  9:  (d)  Temperature,  (e)  u0  +  u.  (f)  v0  +  v. 
Reconstruction  with  the  estimated  systematic  errors  by  Algorithm  3:  (g)  Temperature,  (h)  u0  +  u.  (i)  v0  +  v. 


and  wind  velocity  components,  respectfully,  with  errors  in  two 
travel  times.  Note  the  unrealistic  scale  of  the  fluctuations:  [5, 
50]  °C,  [—33,  8]  m/s,  and  [—17,  16]  m/s.  The  wind  direction, 
depicted  by  streamlines  with  arrows  in  Fig.  6(d),  is  also  signifi¬ 
cantly  distorted.  In  contrast,  using  travel  times  obtained  with 
Algorithm  3,  the  reconstructions  become  relatively  accurate, 
see  Figs.  6(g),  6(h),  and  6(i).  The  magnitudes  of  all  fields  are 
slightly  underestimated,  but  the  main  features  of  the  actual 
fields  are  captured  correctly.  For  example,  in  Fig.  6(g),  the 
locations  and  shapes  of  the  two  eddies  are  correctly  estimated 
with  respect  to  the  actual  temperature  field  shown  in  Fig.  6(a). 
The  similarity  of  the  wind  velocity  components  is  less 
obvious,  but  the  locations  of  “fast”  and  “slow”  eddies  and  their 
orientation  are  captured  with  the  discrepancy  in  the  wind  ve¬ 
locity  magnitude  of  approximately  0.3-0.4m/s. 

V.  EXPERIMENTAL  DATA 

We  now  apply  results  from  the  previous  section  to  ex¬ 
perimental  data  obtained  at  the  BAO.  The  experiment  was 


conducted  on  9  July  2008,  21:32  UTC  at  the  BAO.  Eight 
bend-over  towers  with  three  speakers  and  five  microphones 
were  located  at  a  height  of  9  m,  at  positions  shown  in 


TABLE  V.  Mean  field  reconstructions  of  very  unstable  LES  fields  with 
Algorithm  1,  Algorithm  2,  and  Algorithm  3. 


Case 

Algorithm 

T  (°C) 

u  (m/s) 

v  (m/s) 

Actual 

28.65 

-4.00 

4.00 

Single  error 

Algorithm  1 

26.1  ±0.8 

-2.9  ±0.6 

6.4  ±  0.7 

Algorithm  2 

25.32  ±  0.06 

-2.78  ±  0.05 

6.87  ±  0.05 

Algorithm  3 

28.80  ±  0.06 

-4.10  ±0.04 

3.93  ±  0.05 

Two  errors 

Algorithm  1 

24  ±  1 

—5  ±  1 

5  ±  1 

Algorithm  2 

21.5  ±0.06 

-6.9  ±  0.05 

5.08  ±  0.05 

Algorithm  3 

28.8  ±  0.06 

-4.13  ±0.04 

3.91  ±  0.05 

Large 

Algorithm  1 

26  ±2 

—2  ±  1 

7  ±  1 

errors  ±5  ms 

Algorithm  2 

26.6  ±  0.06 

-2.33  ±  0.05 

5.50  ±  0.05 

Small 

errors  ±  1  ms 

Algorithm  1 

27.0  ±  0.3 

-4.4  ±0.2 

4.2  ±  0.2 

Algorithm  2 

27.51  ±0.06 

-4.22  ±  0.05 

3.90  ±  0.05 
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FIG.  7.  (Color  online)  Reconstruction  the  fields  in  the  BAO  experiment  on  9  July  2008,  21:32  UTC.  Note  significant  improvement  in  reconstructions  (right 
column)  when  Algorithm  3  was  used,  (a)  Temperature  field,  no  estimation  of  the  systematic  errors,  (b)  Temperature  field,  refined  estimation  of  the  systematic 
errors,  (c)  First  horizontal  component  of  the  wind  velocity  vector,  Uq  +  u,  no  estimation  of  the  systematic  errors,  (d)  First  horizontal  component  of  the  wind  ve¬ 
locity  vector,  Uq  +  u,  refined  estimation  of  the  systematic  errors,  (e)  Second  horizontal  component  of  the  wind  velocity  vector,  v0  +  v,  no  estimation  of  the  sys¬ 
tematic  errors,  (f)  Second  horizontal  component  of  the  wind  velocity  vector,  Vq  +  v,  refined  estimation  of  the  systematic  errors. 


Fig.  1(d).  The  systematic  errors  in  the  travel  times  during 
this  experiment  might  have  been  caused  by  errors  in  time 
delays  estimates  and  transducers’  coordinates  that  were  not 
measured  at  the  time  of  the  experiment.  (Such  measurements 
are  costly,  are  time  consuming,  and  are  done  at  the  BAO 

1810  J.  Acoust.  Soc.  Am.,  Vol.  134,  No.  3,  September  2013 


tomography  array  only  from  time  to  time).  Furthermore,  af¬ 
ter  a  previous  experiment,  the  towers  of  the  BAO  array  were 
bent  over  and  then  elevated  again  so  that  the  measured  trans¬ 
ducer’s  coordinates  might  not  exactly  coincide  with  those 
during  the  experiment.  Thus,  the  experiment  on  9  July  2008 
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TABLE  VI.  Mean  field  reconstructions  for  experimental  data.  “Algorithm 
1”  and  “Algorithm  3”  refer  to  the  algorithms  without  systematic  errors  and 
the  refined  reconstruction,  respectively. 


Mean  fields 

T(°C) 

u  (m/s) 

v  (m/s) 

Algorithm  1 

24.0  ±  0.3 

-0.67  ±  0.3 

-0.68  ±  0.3 

Algorithm  3 

25. 14  ±0.09 

-0.71  ±0.07 

-1.11  ±0.07 

presents  a  good  opportunity  to  study  the  reconstruction  of 
temperature  and  wind  velocity  fields  in  the  presence  of  pos¬ 
sible  systematic  errors  in  the  travel  times. 

Table  VI  shows  the  mean  held  reconstructions  by 
Algorithm  1  and  Algorithm  3.  Figures  7(a),  7(c),  and  7(e) 
depict  the  full  reconstructed  temperature  and  wind  velocity 
component  fields,  respectively.  As  one  can  see,  the  recon¬ 
struction  by  the  original  algorithm  looks  unrealistic.  First, 
the  ranges  of  variations  are  [23.0,  25.6]  °C  for  the  tempera¬ 
ture  held  and  [—4.0,  2.5]  m/s  and  ]— 2.60,  1.55]  m/s  for  the 
wind  velocity  components,  which  seem  implausibly  large  for 
the  80  m  by  80  m  tomography  area.  Second,  the  wind  how 
depicted  by  black  lines  in  Fig.  7(a)  is  dramatically  curled, 
which  also  appears  unrealistic.  This  is  reminiscent  of  Figs. 
4(a),  4(c),  and  4(e),  indicating  possible  presence  of  the  sys¬ 
tematic  errors  in  the  measured  travel  times. 

Figure  8  presents  an  analysis  of  systematic  errors.  The 
assessment  of  systematic  errors  assuming  that  they  are  pres¬ 
ent  along  all  paths  (Algorithm  2)  are  shown  in  Fig.  8(a).  A 
more  detailed  analysis  reveals  that  only  six  errors,  corre¬ 
sponding  to  paths  5,  8,  10,  11,  12,  and  15,  have  relatively 
large  magnitudes,  larger  than  0.4  ms.  Following  the  numeri¬ 
cal  case  studies,  the  rehned  reconstruction  by  Algorithm  3 
was  implemented  for  more  accurate  systematic  error  assess¬ 
ments  in  these  paths,  whereas  the  systematic  errors  along 
other  ray  paths  were  set  to  zeros.  The  systematic  errors  for 
this  reconstruction  are  shown  in  Fig.  8(b).  The  mean  held 
reconstruction  by  Algorithm  3  is  shown  in  Table  VI,  while 
the  full  fields  are  depicted  in  Figs.  7(b),  7(d),  and  7(f).  The 
ranges  of  variability  of  the  reconstructed  fields  within  the  to¬ 
mographic  area  have  been  significantly  reduced.  They  are 
[24.9,  25.5]  °C  for  the  temperature  held  and  [—1.65,  0.20] 


m/s  and  [—1.78,  —0.55]  m/s  for  the  wind  velocity  compo¬ 
nents,  which  is  more  realistic.  Also,  the  reconstructed  air 
how  becomes  more  realistic;  see  Fig.  7(d). 

Other  acoustic  tomography  experiments  at  the  BAO 
conducted  at  different  meteo  conditions  (day,  night,  and  dis¬ 
tinct  dates)  were  also  processed.  The  results  obtained  were 
similar  to  those  presented  above,  namely,  the  expected  errors 
decreased  using  Algorithms  2  and  3,  and  the  reconstructed 
helds  had  plausible  ranges  of  variability.  Figure  9  exempli¬ 
fies  a  reconstruction  of  the  helds  on  11  August  2008,  21:31 
UTC.  Note  that  reconstructions  implemented  using 
Algorithm  1  have  unrealistically  wide  ranges  of  variability: 
The  temperature  held  varies  from  21.1  °C  to  26.2  °C  [Fig. 
9(a)],  the  hrst  wind  velocity  component  varies  from  —7  m/s 
to  0.4  m/s  [Fig.  9(c)],  and  the  second  component  lies  in  the 
range  0.8-9. 5  m/s  [Fig.  9(e)].  The  reconstructions  performed 
with  Algorithm  3  are  shown  in  Figs.  9(b),  9(d),  and  9(f), 
respectively.  Note  plausible  ranges  that  became  of  the  same 
order  of  magnitude  as  in  the  LES  numerical  studies.  These 
results  corroborate  those  depicted  in  Fig.  7. 

VI.  SUMMARY  AND  CONCLUSIONS 

Systematic  errors  in  measured  travel  times  can  signifi¬ 
cantly  distort  tomographic  reconstructions.  This  paper  stud¬ 
ied  two  algorithms  aimed  at  assessing  systematic  errors 
present  in  the  data.  Simultaneous  reconstruction  of  the  helds 
and  systematic  errors,  as  considered  previously  in  under¬ 
water  acoustic  tomography,26-28  was  performed.  The  numer¬ 
ical  study  suggests  that  it  can  be  done  successfully  only  if 
the  systematic  errors  are  present  in  a  few,  but  not  all,  of  the 
measured  travel  times.  If  all  travel  times  are  corrupted,  a 
good  reconstruction  is  generally  impossible  unless  the  sys¬ 
tematic  errors  in  travel  times  are  small,  on  the  order  of  1% 
relative  to  travel  times  in  the  considered  experiments.  Based 
on  these  observations,  we  proposed  the  algorithm  for  refin¬ 
ing  reconstructions.  The  initial  systematic  error  estimates 
obtained  with  Algorithm  2  are  used  to  indicate  the  path  indi¬ 
ces  where  the  systematic  errors  are  likely  to  be  present. 
Then,  using  this  information,  systematic  errors,  as  large  as 
6%,  as  considered  in  case  3  of  the  numerical  study,  and 
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FIG.  8.  (Color  online)  Systematic  error  analysis  for  the  BAO  experiment  on  9  July  2008,  21:32  UTC.  (a)  Systematic  errors  estimated  by  Algorithm  2,  assum¬ 
ing  systematic  errors  in  all  paths,  (b)  Systematic  errors  estimated  by  the  refined  Algorithm  3,  assuming  systematic  errors  in  paths  5,  8,  10,  11,  12,  and  15. 
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FIG.  9.  (Color  online)  Reconstruction  the  fields  in  the  BAO  experiment  on  11  August  2008,  21:31  UTC.  Note  significant  improvement  in  reconstructions 
(right  column)  when  Algorithm  3  was  used,  (a)  Temperature  field,  no  estimation  of  the  systematic  errors,  (b)  Temperature  field,  refined  estimation  of  the  sys¬ 
tematic  errors,  (c)  First  horizontal  component  of  the  wind  velocity  vector,  u0  +  u,  no  estimation  of  the  systematic  errors,  (d)  First  horizontal  component  of  the 
wind  velocity  vector,  u0  +  u,  refined  estimation  of  the  systematic  errors,  (e)  Second  horizontal  component  of  the  wind  velocity  vector,  v0  +  v,  no  estimation  of 
the  systematic  errors,  (f)  Second  horizontal  component  of  the  wind  velocity  vector,  Vq  +  v,  refined  estimation  of  the  systematic  errors. 


mean  fields  can  be  estimated  more  accurately  by  setting  the 
remaining  errors  to  zero  (Algorithm  3). 

The  developed  algorithms  were  studied  numerically  on 
several  LES  data  sets  corresponding  to  very  unstable,  unsta¬ 
ble,  and  neutral  atmosphere,  and  then  applied  to  experimental 
data  acquired  at  different  conditions  (day  and  night),  although 
typical  results  and  observed  trends  pertinent  to  two  numerical 
and  outdoor  data  sets  were  reported  in  the  paper.  The  results 
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obtained  showed  that  expected  errors  decrease  remarkably  in 
the  numerical  and  outdoor  experiments  and  the  reconstruc¬ 
tions  became  more  plausible,  in  the  meaningful  physical 
ranges  for  a  given  stability  regime  in  the  atmosphere. 

As  the  suggested  algorithms  for  assessment  and  elimina¬ 
tion  of  systematic  errors  are  general,  they  can  be  applied  for 
other  outdoor31  and  indoor12  acoustic  tomography  arrays. 
However,  particular  results  of  this  approach  might  depend 
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on  the  geometry  of  the  transducers’  locations  in  a  tomogra¬ 
phy  array.  Therefore,  further  studies  are  necessary  to  assess 
effects  for  different  arrays,  their  geometry  and  sizes,  and 
fraction  of  the  error-free  measurements  for  successful  appli¬ 
cation  of  the  considered  algorithms. 
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